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A recently developed model for the QCD analytic invariant charge is compared with 
quenched lattice simulation data on the static quark-antiquark potential. By employing 
this strong running coupling one is able to obtain the confining quark-antiquark potential 
in the framework of the one-gluon exchange model. To achieve this objective a technique 
for evaluating the integrals of a required form is developed. Special attention is paid here 
to removing the divergences encountered the calculations. All this enables one to examine 
the asymptotic behavior of the potential at both small and large distances with high 
accuracy. An explicit expression for the quark-antiquark potential, which interpolates 
between these asymptotics, and satisfies the concavity condition, is proposed. The derived 
potential coincides with the perturbative results at small distances, and it is in a good 
agreement with the lattice data in the nonperturbative physically— relevant region. An 
estimation of the parameter Aq CD is obtained for the case of pure gluodynamics. It is 
found to be consistent with all the previous estimations of Aqcd in the framework of 
approach in hand. 
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1. Introduction 

The description of hadron dynamics at large distances remains a crucial challenge of 
elementary particle physics for a long time. On the one hand, the asymptotic free- 
dom of Quantum Chromodynamics (QCD) allows one to apply perturbation theory 
to study some "short-range" phenomena, for instance, the lowest-lying bound states 
of heavy quark-antiquark systems. On the other hand, for the consistent description 
of many phenomena related to the "long-range" dynamics (such as confinement of 
quarks, structure of the QCD vacuum, etc.) something more than the usual pertur- 
bative approach has to be involved. 
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Theoretical analysis of the strong interaction basically relies upon the renor- 
malization group (RG) method. Usually, in order to describe the hadron dynamics 
in the asymptotical ultraviolet (UV) region, one applies the RG method together 
with perturbative calculations. However, this leads to unphysical singularities of 
the outcoming solutions to RG equations, that contradicts the general principles 
of local Quantum Field Theory (QFT). An effective way to overcome such diffi- 
culties consists in invoking the analyticity requirement. This pres cription became 
the underlying idea of the so-called analytic approach to QFTjJ*™ which was lately 
extended to QCD by Shirkov and Solovtsov. L " u - ±l 

In the framework of the analytic approach to Q CD a new model for the strong 
running coupling has recently been developed!"^ Its basic idea is to impose the 
analyticity requirement on the RG (3 function perturbative expansion for restoring 
its correct analytic properties (see Sec. 2 for the details). The analytic invariant 
charge possesses a number of profitable properties, in particular, it contains no 
unphysical singularities (see also Refs.0IHlElfbr the description of its traits). It is 
worth noting here that this analytic running coupling incorporates the perturbative 
and intrinsically nonperturbative features of the strong interaction, and enables one 
to describe a wide range of hadron processes. 

A crucial insight into the nonperturbative aspects of the strong interaction at 
large distances can be provided by lattice simulations. This may be, for instance, 
the stati c qu ark-antiquark (QQ) potential calculated up to the significantly large 
distances^ (r > 1 fm), or the investigation of the topological structure of the QCD 
vacuum \^ Certainly, a decisive test of any model for the strong interaction is its 
comparison with the lattice results. 

The objective of this paper is to study the asymptotic behavior of the static 
quark-antiquark potential constructed by making use of the analytic invariant 
charge. It is also of a primary interest to collate the derived QQ potential with 
the relevant lattice simulation data. For verification the consistency of the results 
it is worth to evaluate the corresponding parameter A QCD and to compare it with 
its previous estimations obtained in the framework of approach in hand. 

The layout of the paper is as follows. In Sec. 2 the model for the QCD analytic 
invariant charge is briefly described and its basic features are outlined. In Sec. 3 the 
asymptotic behavior of the quark-antiquark potential, constructed by making use of 
this strong running coupling, is examined. A technique for treating the divergent in- 
tegrals encountered is developed here. The explicit expression for the QQ potential, 
which possesses the obtained asymptotics, and satisfies the concavity condition, is 
derived. In Sec. 4 these results are applied to study of the relevant lattice simulation 
data. The constructed quark-antiquark potential agrees fairly well with the lattice 
data in the nonperturbative physically-relevant region. At the same time, it coin- 
cides with the perturbative QQ potential at small distances. The fitted parameter 
Aqcd is found to be in a good agreement with all its previous estimations, implying 
the consistency of obtained results. In the Conclusion (Sec. 5) the achieved goals 
are formulated in a compact way. 
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2. The QCD Analytic Invariant Charge 

As it has been mentioned in the Introduction, the perturbative approximation of the 
(3 function in the RG equation for the strong running coupling a(/j, 2 ) = g 2 (n 2 ) / (Air) 

leads to unphysical singularities of the outcoming solutions, that contradicts the 
general principles of local Quantum Filed Theory. A plausible way to overcome 
such difficulties consists in imposing the analyticity requirement on the perturba- 
tive expansion of the RG [3 function for restoring its correct analytic properties. This 
prescription is a distinctive feature of the recently developed model for the QCD an- 
alytic invariant charge02l At the one-loop level the corresponding renormalization 
group equation can be solved explicitly: 

Po zmz A z 
while at the higher loop levels only the integral representation for the analytic 
invariant charge has been obtained (see Refs. H3 EJl- Figure presents the analytic 
running coupling 5 an (<Z 2 ) = a a n(<7 2 )/3o/(47r) at different loop levels. 
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Fig. 1. The QCD analytic invariant charge at different loop levels, z = q 2 /A 2 . 

The developed model for the QCD analytic invariant charge possesses a num- 
ber of appealing traits (see Refs. H3 E|). Namely, it has no unphysical singularities 
at any loop level; it contains no adjustable parameters; it incorporates ultraviolet 
asymptotic freedom with infrared (IR) enhancement in a single expression; it has 
universal behavior both in UV and IR regions at any loop level; it possesses a good 
higher loop and scheme stability. This model also enables one to describe various 
strong interaction processes both of perturbative and intrinsically nonperturbative 
nature (see papers 171 1121 1181 and references therein). It is of a primary importance to 
mention here that the analytic invariant charge (J2J has recently been rediscovered 
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when studying the conformal inversion symmetry related to the size distribution of 
instantonsM^I In particular, Eq. was proved to reproduce explicitly this kind of 
symmetry. In turn, the latter is i n a good agreement with the relevant lattice data 
by the UKQCD Collaboration^] ( see Refs. [Hand for the details). 



3. The Static Quark Antiquark Potential 

Let us proceed now to the construction of the static quark-antiquark potential. In 
the framework of the one-gluon exchange model it is related to the strong running 
coupling a(q 2 ) by the 3-dimensional Fourier transformation 



Strictly speaking, this definition of the potential is justified for small distances 
(r < 0.1 fm) only. For example, the lowest-lying bo und states of heavy quarks can 
be described by employing the perturbative a QCD|l2l However, at large distances 
(r > 0.5 fm), which play the crucial role in hadron spectroscopy, the perturbative 
approach blows up due to unphysical singularities (such as the Landau pole) of the 
strong running coupling. Nevertheless, the model J2Jl, being complemented with a 
certain insight into the nonperturbative behavior of the QCD invariant charge, has 
proved to be successful for description both heavy-quark and light -quark systems 
(see, e.g., reviews ITffl II 71 ITH1 and references therein for the details). 

In this paper, for the construction of the static potential of the quark-antiquark 
interaction, we shall use the analytic invariant charge 10) ■ After integration over 
angular variables, Eq. © in this case takes the form 



v y) = ~w — / ~T~i — 2 — d P> 4 
3/3 r J p 2 lap 2 pR 



where p — qr , R = r/r , and r is a reference scale of the dimension of length, 
which will be specified below. This integral diverges at the lower limit, that is a 
common feature of the models of such kind (see, e.g., Ref. We shall treat this 
divergence by employing the analytical rcgularization (see Eq. QDj) ). 

The analysis of Eq. performed in Ref. 5, elucidated only the leading asymp- 
totic behavior of V(r) at large distances. However, it is still desirable to examine 
at a more precise level the asymptotics of the quark-antiquark potential Q in the 
both, ultraviolet and infrared regions. This objective can be achieved in the follow- 
ing way. First of all, let us introduce a dimensionless variable Q = pR and rewrite 



a The leading short-distance nonperturbative effect due to the gluon condensate has also been 
taken into account in Ref. 1151 
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Eq. gj: 



v(r) - 16 1 1 f(i-^\ dQ 



(5) 



Then, formally expanding the denominator of the integrand, 13 one can find the 
asymptotic behavior of V(r) at both small and large distances: 



V(r) 



16 1 



1 n ° 1 

R ^ (lni?)" +1 

n=0 y ' 



sinQ 



(InQ)™ dQ 



o 

sin Q 



(6) 



The values of no and mo will be specified below. 

For evaluation of the expansion coefficients in Eq. |0 it is worth considering an 
auxiliary integral of the form: 



_ , r(i + |1 

Q* sinQdQ = V^2* — ^ |4 



(7) 



Differentiating this equation n times with respect to variable t we obtain 



Q t (lnQ) n smQdQ = v(n,t), 



(8) 



where 



d n 

v(n, t) — \TH — — 



r(i 



(9) 



In Eq. (0 the parameter t takes the values < |Re (t + 1)| < 1 (see, e.g., Ref. EOJ. 
However, in order to determine the expansion coefficients of the second sum in 
Eq. ©, one has to go to the point t = —3, which is outside of this range. Neverthe- 
less, it can be done by making use of the analytic continuation of the left-hand side 
of Eq. JJJ • This continuation is unique and it is defined obviously by the right-hand 
side of Eq. JJJ) all over the complex i-plane except for the points t = —2N, with N 
being a natural number. Fortunately, we are not dealing with these values of the 
parameter t, and we can put 



sin Q 



(lnQ)™ dQ = v(n,t) 



s^2N (TV = 1,2, 3,. 



(10) 



b A similar method has also been used in Ref. 1191 
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It should be noted here that this analytical continuation plays the role of regular- 
ization of the expansion coefficients in Eq. JfjJ). It is also convenient to introduce the 
notations 



■ v(n, t) 



v(m, t) 



(11) 



A few first coefficients 1)110 have a quite simple form, namely u a = 1, Ux = —7, 
u 2 = 7 2 + tt 2 /12; and uj = 1/2, Ui = 3/4 - 7 /2, oj 2 = j 2 /2 - 3 7 /2 + tt 2 /24 + 7/4, 
where 7 ~ 0.57721... denotes the Euler's constant (see Ref. 0. It is worth noting 
here that the leading coefficients uq, ojo, and oji have also been calculated in Ref.El 
but by making use of another technique. 

Thus, the static quark-antiquark potential © can be represented now in the 
following form: 



V(r) - I 
1 ' 3/3 r 



"0 



^ (lni?)'^ 1 

n=0 v ' 



R 



mo 



^ (ln J R)™+ 1 

m=0 v ' 



R=-. (12) 



At small distances the potential l|12|) possesses the standard behavior, determined 
by the asymptotic freedom 

.?r 1 1 



V(r) 



0. 



3/3 r i? lni?' 
At the same time, it proves to be rising at large distances 

R 



rri s 87T 1 

V(r) = — 



(13) 



(14) 



3/3 r 2 lni?' 

implying the confinement of quarks. It is of a particular interest t o m ention here 
that a similar rising behavior of the QQ potential has been proposed^ a long time 
ago proceeding from the phenomenological assumptions. 

Equation lfT2^l describes the behavior of the static quark-antiquark potential (@J 
at small and large distances. However, its straightforward extrapolation to all dis- 
tances encounters poles of different orders at the point R = 1, which apparently is 
an artifact of the expansion For the practical purposes it would be undoubtedly 
useful to derive an explicit expression for the QQ potential applicable for < r < 00. 

In order to construct an explicit interpolating formula for the quark-antiquark 
potential (@J we shall employ here the following method. Let us modify the expan- 
sion H12JI in a "minimal" way, by adding terms which only subtract the singularities 
at the point R — 1 and do not contribute to the derived asymptotics. This leads to 
the following expression for the static quark-antiquark potential 



V(r) 



1 



RQnRY 



Reg 



E 



R 


Reg} 







(15) 



where R = r/r , and at the leading orders [1/(R lni?)]R Cg = 1/ (R In R) — 1/P, 
[l/(i?(lni?) 2 )] Rcg = 1/(RQjiR) 2 ) - 1/P 2 , and [R/lnR} Rcg = R/JnR- 1/P, 
[i?/(lni?) 2 ] Rcg = R/{\nR) 2 - 2/P - 1/P 2 , P = R-1 (see also Ref.0). 
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Fig. 2. The quark— antiquark potential 1151 in dimensionless units at different levels of approx- 
imation: no = mo = (dashed curve), no = mo = 4 (solid curve), no = mo = 5 (A), and 
no = mo = 10 (■). The one-loop perturbative result is shown by dotted curve. 



The numerical analysis of Eq. (|15|l revealed that for practical purposes it is 
enough to retain the first five expansion terms (no = too = 4) therein. It turns 
out that the potential Ijl5(l itself and the corresponding estimation of the value of 
parameter A QCD are not affected by higher-order contributions. In particular, the 
curves (|15|) for no — too — 5 and for no — too = 10 are practically indistinguishable 
of the curve corresponding to no = too = 4 over the whole region < r < oo (see 
Figure E). And the higher-order estimations of the parameter A QCD vary within 
0.5 % of the value obtained at the no — too = 4 (see Sec. 4 further). 

Thus, we arrive at the following explicit expression for the static quark-antiquark 
potential: 



V(r) = Vo + 



0.577 




1.156 

(lni?) 3 

3.185 



4.021 
" (lni?) 4 

17.844 



15.018 

(lni?) 5 



(lni?) 4 (lni?) 5 
84.143 32.861 



(1-i?) 2 (1-i?) 3 (1-i?) 4 (l-i?) £ 



(16) 



where i? = r/ro- In order to reproduce the correct short-distance behavior of the 
QQ potential, the dimensional parameter ro in this equation has to be identified with 
Ams by the relation r^ 1 = Aexp(7) (see, e.g., Refs. 1221 1231 124fl . It is straightforward 
to verify that the potential (|16fl satisfies also the concavity condition 

d 2 V{r) 



dV(r) 



>0, 



< 0, 



dr dr 2 
which is a general property of the gauge theories (see Ref. 



(17) 



| for the details). 
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4. Discussion 

In this section we are going to apply the obtained results to the study of recent 
quenched lattice simulation datcM^on the static quark-antiquark potential. In par- 
ticular, the value of the parameter A QCD will be estimated by comparing the derived 
expression for the QQ potential 116f> with these lattice data. 




Fig. 3. Comparison of the quark-antiquark potential V(r) defined by Eq. 1161 (solid curve) with 
the quenched lattice simulation datalSl (») . The values of the parameters are: A = 670 MeV, rif = 0, 
Vo = —3.164 GeV. The dashed curve corresponds to the relevant one-loop perturbative result. 



Thus, a fit of the quark-antiquark potential 1|16|) to the quenched lattice sim- 
ulation datcP^ has been performed with the use of the least square method. The 
varied parameters in Eq. Ijl6(l were r$ and the additive self-energy constant Vo- 
The result of the fit is presented in Figure This figure shows that in the non- 
perturbative physically-relevant range 0.3 fm < r < 1.2 fm, in which the average 
quark separations \J (r 2 ) for quarkonia sitsj^the QQ potential {EH) reproduces the 
lattice datcP^I fairly well. At the same time, in the region r < 0.05 fm the derived 
potential coincides with the perturbative result (dashed curve in Fig. EJ . The dif- 
ference between the lattice data and the expression l|16[) in the intermediate range 
0.05 fm < r < 0.3 fm may be explained by the presence of additional nonpertur- 
bative contributions at these distances (see, e.g., Refs. 0§1 130fl . But the detailed 
investigation of this matter is beyond the scope of the paper. 

Thus, the estimation of the parameter A QCD in the course of this comparison 
gives A(" f=0 ) = (670 ± 8) MeV (this value corresponds to the one-loop level with 
n f = active quarks). The uncertainty here has been calculated by making use 
of the "3cr-criterion" (see Ref. 13 It . It is worth noting that the fit with the use of 
the maximum likelihood method results in even better reproduction of the lattice 
data and gives a similar value A(™ f=0 ) ~ 645 MeV. Further, in order to collate the 

c The reliability of the perturbative quark-antiquark potential at small distances has been discussed 
in Refs. l23ll27ll28l 
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obtained estimation with the earlier ones, one has to continue it to the region of 
three active quarks. It was performed by employing the matching procedure, and 
gives d the value A(™ f=3 ) = (590 ± 10) MeV. The la tter is in a good agreement with 
all the previous estimations of this parameteJSEH [A(™ f=3 ) = (550 ± 50) MeV] in the 
framework of the approach in hand. 

5. Conclusion 

In the paper the QCD analytic invariant charge has been applied to study of the 
quenched lattice simulation data on the static quark-antiquark potential. This 
strong running coupling enables one to obtain explicitly the confining quark- 
antiquark potential in the framework of the one-gluon exchange model. A technique 
for evaluation the integrals of a specific form, developed in this paper, allows one 
to examine the asymptotic behavior of the derived potential at small and large dis- 
tances with high accuracy. An explicit formula for the quark-antiquark potential, 
which interpolates between these asymptotics and satisfies the concavity condition, 
is obtained. The derived potential has the standard form determined by asymptotic 
freedom at small distances. At the same time, in the nonperturbative physically- 
relevant region 0.3 fm < r < 1.2 fm this potential agrees fairly well with the lattice 
simulation data. The value of the parameter A QCD is estimated for the case of pure 
gluodynamics. It is found to be consistent with all the previous estimations of A QCD 
in the framework of the approach developed. 
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